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BASE  MOTION  ISOLATION  OF  A 
TWO-AXIS  BEAM  DEFLECTOR 
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ABSTRACT 

■  The  base  motion  isolation  characteristics 
of  an  optical  beam  deflector  with  a  two-axis 
gimbal  support  are  investigated.  The  particular 
configuration  of  the  beam  deflector  system  con¬ 
sidered  utilizes  rate  integrating  gyros  to 
provide  a  rate  inner  control  loop  with  an 
optical  tracker  providing  an  outer  position 
control  loop.  Base  motion  isolation  is  par¬ 
tially  achieved  by  the  gimbal  support  assumed 
to  be  frictionless.  However,  input  beam  motion 
relative  to  the  optical  surface  of  the  beam 
deflector  is  shown  to  produce  large  excursions 
of  the  output  beam.  As  a  result,  a  feed-forward 
control  loop  is  added  to  measure  and  compensate 
for  input  beam  motion,  i.e.,  base  motion.  Every 
effort  is  directed  to  use  optimistic  estimates 
of  beam  deflector  performance  in  order  to 
develop  a  measure  of  the  smallest  pointing 
error  possible  with  such  a  two-axis  beam 
deflector  configuration.  The  base  motion 
spectrum  used  is  intended  to'  approximate  that 
of  a  large  aircraft  in  normal  flight. 


I.  INTRODUCTION 

Recently  applications-  have  arisen  which 
involve  the  airborne  transmission  and/or 
reception  of  information  with  optical  wave¬ 
lengths.  Many  of  these  applications  require  a 
precise  optical  pointing  system  to  transmit  or 
collect  the  optical  energy.  Several  configura¬ 
tions  of  optical  surfaces  are  conceivable  for 
such  a  purpose.  One  of  the  simplest  concepts 
to  conceive  is  that  of  the  optical  beam 
deflector  as  depicted  in  figure  1.  In  order 
to  provide  the  capability  of  directing  the 
optical  beam  in  any  direction  independent  of 
aircraft  altitude,  it  is  necessary  to  mount 
the  beam  deflector  in  gimbals.  The  particular 
gimbal  configuration  depicted  in  Figure  2 
allows  for  motion  of  the  beam  deflector  about 
two  axes.  This  is  the  minimum  number  of 
gimbals  needed  to  allow  motion  in  both  azimuth 
and  elevation. 


Figure  2.  Two-Axis  Beam  Deflector 
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Depending  upon  the  specific  application  of 
the  optical  Information,  the  beam  deflector  may 
be  used  either  to  transmit  optical  energy  from 
the  aircraft  to  a  target  or  to  collect  optical 
energy  from  an  external  target  source  into  the 
aircraft.  For  some  applications  two-way  trans¬ 
mission  of  optical  energy  might  be  required. 

For  purposes  of  clarity  only,  the  source  of 
optical  energy  will  be  considered  to  be  aboard 
the  aircraft  and  the  beam  deflector  will  have 
the  purpose  of  directing  this  energy  upon  an 
external  target.  Thus  the  input  beam  will  be 
considered  to  be  the  optical  energy  transmitted 
from  the  airborne  source  to  the  beam  deflector 
while  the  output  beam  will  be  the  optical  energy 
path  from  the  beam  deflector  to  the  external 
target.  For  applications  in  which  the  role  of 
source  and  target  are  interchanged  the  follow¬ 
ing  investigation  applies  fully  with  only  a 
change  in  nomenclature. 

The  investigation  of  base  motion  isolation 
capability  of  the  particular  beam  deflector 
configuration  considered  herein  is  intended  to: 

a)  determine  the  smallest  pointing  error 
which  could  conceivably  be  obtained  from  such  a 
configured  system  operating  within  an  aircraft 
with  the  base  motion- spectrum  assumed,  and 

b)  serve  as  a  guide  to  considering  the 
effects  of  base  motion  upon  the  performance  of 
och-v  airborne  optical -pointing  systems. 


II.  3EAM  DEFLECTOR  COHTROL 

In  order  to  discuss  quantitatively  the 
control  and  performance  of  the  two-axis  beam 
deflector,  it  is  necessary  to  introduce  several 
coordinate  systems.  These  are  depicted  in 
Figure  3  as: 

I  -  fixed  in  inertial  space 

V  -  fixed  to  the  vehicle  with  coordinate  axes 

i,  j,  h 

0  -  fixed  to  the  outer  or  azimuth  gimbal  with 
coordinate  axes  n,  e,  k 

M  -  fixed  to  the  beam  deflector  with  coordinate 
axes  r,  e,  d 

L  -  line-of-sight  coordinates  with  coordinate 
axes  1,  e,  p. 


Figure  3.  Coordinate  Axes 


The  purpose  of  the  beam  deflector  control 
system  is  to  direct  the  output  beam  along  the 
desired  direction.  This  is  accomplished  by 
controlling  the  elevation  angle  e  and  the 
azimuth  angle  n  where  n  is  measured  between 
the  base  and  the  outer  gimbal  and  c  is  measured 
between  the  outer  gimbal  and  the  beam  deflector 
surface  as  shown  in  Figure  A. 


-  2  - 


Figure  4.  Controlled  Angles 


The  beam  deflector  control  can  be  accom¬ 
plished  by  two  control  loops,  elevation  and 
azimuth.  These  basic  features  of  these  loops 
are  block  diagrammed  in  Figures  5a  and  5b. 


For  the  particular  configuration  of  beam 
deflector  considered  herein,  it  is  assumed  that 
the  input  beam  source  is  mounted  on  the  outer 
gimbal  coincident  with  that  gimbal's  axis  of 
rotational  freedom.  For  such  a  configuration, 
an  angular  velocity  of  the  outer  gimbal  and 
hence  the  input  beam,  about  the  elevation  axis, 
“0  >  produces  a  corresponding  angular  velocity 
e 

of  opposite  sign  of  the  output  beam  direction. 
This  is  accounted  for  in  Figure  5a  by  the  dis¬ 
turbance  angular  velocity,  -uQ  . 

e 

A  rotation  of  the  mirror  surface  with 
respect  to  the  input  direction  about  the  e  axis 
is  easily  seen  to  produce  a  rotation  of  twice 
that  magnitude  of  the  output  beam  direction. 


The  angular  velocity  of  the  output  direction 
about  the  e  axis  is  thus 

UL  "  2“m  '  “0  (1) 

e  e  e 

Since  the  direction  to  be  controlled  is 
that  of  the  output  beam,  it  is  desirable  to 
measure  the  motion  of  this  output  direction 
with  a  high  bandwidth  sensor.  Unfortunately, 
the  output  direction  is  only  a  direction  in 
space;  it  does  not  correspond  to  any  direction 
in  the  physical  hardware.  Thus  there  is  no  way 
of  mounting  a  sensor  on  this  output  beam  direc¬ 
tion  to  measure  its  motion.’*'  Instead  motion 
sensors  can  only  be  attached  to  the  gimbals  and 
the  vehicle  itself.  The  elevation  control  loop 
of  Figure  5a  postulates  the  use  of  a  rate 
integrating  gyroscope  mounted  on  the  inner 
gimbal  with  its  input  direction  along  the  e  axis. 
The  azimuth  control  loop  of  Figure  5b  utilizes  a 
rate  integrating  gyroscope  mounted  on  the  outer 
gimbal  with  its  input  direction  along  the.  k  axis. 
The  rate  integrating  gyros  provide  a  relatively 
high  bandwidth  rate-feedback  loop  about  the 
elevation  and  azimuth  gimbals. 

Inasmuch  as  the  measurement  of  primary 
Interest  is  the  orientation  of  the  output  berm 
with  respect  to  the  target,  an  optical  tracker 
is  incorporated  in  the  control  system.  The 
tracker  is  oriented  so  that  the  center  of  its 
field  of  view  always  is  coincident  with  the 
output  beam  direction.  Thus  the  tracker 
measures  angles  between  the  output  beam  direc¬ 
tion  and  the  desired  target  direction.  The 
tracker  thus  provides  an  outer  position-feedback 
loop. 

The  beam  deflector  is  positioned  by  gimbal 
torque  motors.  The  conmands  are  input  to  the 
torque  generator  of  the  rate  integrating  gyro¬ 
scope.  The  output  angle  of  the  gyroscope  pro¬ 
duces  the  electrical  input  to  the  gimbal 
torque  motor. 

Figure  5  depicts  disturbance  torques  T^ 

e 

and  Td  upon  the  elevation  and  azimuth  control 

“a 

loops  respectively.  These  torques  are  dynamic 
reaction  torques  arising  from  the  rotation  of 


*  Other  mechanizations  of  the  two-axes  beam 
deflector  are  possible  in  which  a  slaved  plat¬ 
form  is  driven  by  the  angular  velocity  of  the 
mirror  with  a  two-to-one  gear  ratio.  The 
slaved  platform  body  axes  then  correspond  with 
those  of  the  line-of-sight  frame.  Measure¬ 
ments  of  the  motion  of  the  slaved  platform  made 
by  sensors  which  can  be  physically  attached  to 
it,  can  thus  be  used  as  measurements  of  the 
motion  of  the  output  beam  direction.  Such  a 
system  has  been  described  elsewhere  ^  and  is 
not  considered  in  this  investigation. 


Figure  5a.  Basic  Elevation  Loop 
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Figure  5b.  Basic  Azimuth  Loop 


the  gimbal  Frames  with  respect  to  inertial 
space.  These  dynamic  torques  have  been  pre¬ 
sented  elsewhere  ^  .  Appendix  A  contains  a 
derivation  of  the  equations  of  motion  for  both 
the  elevation  and  azimuth  gimbals  Including 
the  dynamic  reaction  torques. 


III.  BASIC  ELEVATION  CONTROL  LOOP  DESIGN  AND 
EVALUATION 

Figure  5a  presents  a  functional  block  dia¬ 
gram  of  a  basic  elevation  control  loop.  In 
order  to  evaluate  the  pointing  performance  of 
such  a  control  loop  the  following  representative 
parameter  values  were  used  s 


■  -  132000 


<£•  +  «  (fy  +  1) 

G(s)  .  - £Z - 

(s  +  1)  <^§2  +  D 

'  In  order  to  track  a  target  with  a  constant 
line-of-sight  velocity  without  a  steady  state 
tracking  error,  a  type  IX  control  loop  is 
desired.  This  is  obtained  by  choosing  the 
tracker  dynamics  as 


The  overall  loop  gain  can  be  adjusted  by 
setting  the  tracker  gain,  Kc.  In  order  to 

provide  a  basis  for  this  choice,  a  root  locus 
of  the  control  loop  with  parameter  values  as 
given  above  is  presented  as  Figure  6. 

To  minimize  the  effect  of  the  disturbances 
produced  by  the  dynamic  reaction  torque  and  by 
the  base  motion,  a  large  value  for  loop  gain  is 
desired.  However,  as  evidenced  by  Figure  6, 
too  large  a  value  for  loop  gain  will  produce  an 
underdamped  response  which  would  lead  to  unaccept¬ 
able  settling  following  commanded  angular  slews. 
Thus  a  compromise  value  of  loop  gain  is  needed. 

The  closed  loop  poles  indicated  on  Figure  6 
result  from  the  choice  of  tracker  gain,  K  ,  of 
500.  " 

The  basic  control  loop  design  is  thus 
completed. 

To  evaluate  the  effectiveness  of  the  basic 
control  loop  design  in  rejecting  undesired  dis¬ 
turbances,  one  must  consider  the  two  sources  of 
disturbance:  the  dynamic  torque,  T^  ,  and  the 

e 

input  beam  motion,  Oq  .  Both  of  these  distur- 
e 

bonces  can  be  written  in  terms  of  the  base  motion 
disturbances  and  the  required  tracking  rate.  In 
functional  form 

Td  “  \  (UV  '  “V  V  %  V'  V  WL> 

e  eijpipjp 


figure  6.  Elevation  Root  Locus 
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u„  is  the  component  of  base  angular 
Vi 

velocity  along  the  i  direction. 

u„  is  the  component  of  base  angular 

J 

velocity  a^ong  the  j  direction. 


is  the  component  of  line  of  sight 
P 

angular  velocity  along  the  p  direction. 


Tne  angular  velocities  of  the  base  and 
line  of  sight  are  predictable  only  as  random 
variables .  Thus  these  may  be  described 
statistically  in  terms  of  their  power  spectral 
density  (PSD)  functions.  The  result  of  the 
analysis  will  be  the  RMS  pointing  error  angle, 
6  .  Since  6  is  the  output  of  a  linear 

CRMS  e 


system  having  several  random  disturbance  Inputs 
assumed  to  be  statistically  independent, 


2 

♦  («) 

Xi 


where  — (w)  is  the  performance  function 
X1 

relating  the  1th  independent  disturbance, 
to  the  pointing  error,  0g 


Tne  RMS  pointing  error  is  then  obtained  as 


0 

CRMS 

In  order  to  carry  out  the  above  evaluation 
procedure,  it  is  necessary  to: 

a.  Write  the  disturbance  inputs  in  terms 
of  the  basic  variables  of  base  motion  and 
pointing  rates. 

b.  Determine  the  PSD  of  each  of  these 
disturbance  inputs  assuming  independence. 

c.  Determine  the  performance  function 
relating  each  of  these  disturbance  inputs  to 
the  pointing  error  angle,  6fi. 

Stepn  a  and  b  above  are  carried  out  in 
Appendix's.  Step  c  is  readily  accomplished 
noting  fi:om  Figure  5a  that 
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The  PSD  chosen  to  represent  aircraft 
angular  velocity  about  each  axis  is  given  as 

„  (0.1)  (1400*) _  (rad/sec)2 


“v  (u2  +  1.52)(io2  +  14002) 


hz 


(2) 


This  base  motion  PSD  is  plotted  as  Figure  7. 

e 

Using  the  procedure  outlined  above  the' 
PSD  of  resulting  pointing  error  for  the  basic 
elevation  control  loop  is  determined.  This 
PSD  is  plotted  as  Figure  8.  The  RMS  pointing 
‘error  angle  is  evaluated  to  be 


9  -  942.8  urad 

eRMS 


Figure  7.  Base  Motion  PSD 


Figure  8.  Pointing  Error  PSD 


.  As  evidenced  by  Figure'S,  the  pointing 
error  arises  due  ,to  both  dynamic  torque  distur¬ 
bance  and  also  due  to  the  angular  velocity  of 
the  input  beam.  .  produced  by  the  base  ( 

motion.  In  'order  to  evaluate  the  relative 
importance  of  each  of  these  disturbances,  the 
pointing  error  resulting  from  the  dynamic  torque 
disturbance  alone  'is  evaluated  using  the  above 
procedures.  The  RMS  pointing  brror  caused  by 
the  dynamic  torque  alone  is  found  to  be  ; 

9  (byt  dynamic  torque)  «•  1.8  prad 
RMS  • 

.  From  the  above,  it  is  readily  concluded 
that  the  overwhelming  cause  of  the  pointing  ! 
error  determined  above  is  base  motion  distur¬ 
bance.  For  this  reason,  the -dynamic  torque 
disturbance  effects  will  not  be  considered 
further. 

For  most1  precision  pointing  applications, 
the  RMS  pointing  error  angle  of  about  1  mijlli- 
radlan  determined;  above  is i not  acceptable.  In 
order  to  reduce  the  pointing  error,  base  motion 
compensation  is  evaluated.  . 


IV.  FEED-FORWARD  COMPENSATION  FOR  BASE  MOTION 

■  i  ■ 

In  order  to  reduce  the  pointing  error 
caused  by  phe  motion  of  the  input  beam,  i.e., 
by  base  motion,  a  feed-foryard  compensation 
scheme  is  considered.  This  scheme  proposes  to 
add  a  compensation  signal,  Xc>  as  a  command 

torque  to  the  .rate  integrating  gyro  in  the, 
elevation  control  scheme  innpr  loop.  The  com¬ 
pensation  command  must  be  chosen  sue!)  that  the 
disturbance  of  lnpht  beam  motion,  -uQ  ,  shown  | 

,  e 

in  Figure  5a  is  effectively  cancelled. 

i 

Consider  the  portion  of  the  elevation 
control  loop  reproduced  from  Figure  5a  as 
Figure  9.  1  !  , 


In  order  to  determine  the  proper  cqmpensa- 
tion  signal,  );J,,  the  performance  function  of 

the  inner  lqop  must  be  Investigated. 

!  "  . 
l(sing  the  dynamics  postulated  previously1, 
the  open  loop  performance  function  of  this 
inner  loop  is 


PF0L<8> 


132000  (^  +  1)  <~  +  1) 
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The  resultant  closed  loop  performance  function 
for  tlus  inner  loop  is 
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1‘ 
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Althoug!,  r.he  above  shows  thst  complete 

compensation  ia  theoretically  possible,  one 

must  consider  the  practical  feasibility  of 

determining  the  locations  of  poles  and  zeroes  - 

of  an  actual  hardware  inner  loop  in  order  to 

construct  G  (s) .  In  reality,  'tne  could  not 
comp 

determine  the  locations  of  singularities  to 
five  decimal  places  as  assumed  ,.oove .  The 
approximate  frequency  response  of -the  actual 
Inner  loop  could  be  determined  over  the 
frequencies  of  Interest  by  exhaustive  testing 
using  known  ir.pt  ts.  The  approximate  locations 
of  the  singularities  could  then  be  determined 
analytically  to  match  the  measured  frequency 
response  in  a  satisfactory  manner.  The 
proximity  of  the  singularities  determined  in 
such  a  manner  to  those  of  the  actual  hardware 
control  loop  would  depend  upon  many  factors 
including  the  accuracy  with  which  the  experi¬ 
mental  measurements  are  obtained  and  the  uktlx 
of  the  investigator  in  choosing  an  analytical 
function  with  the  correct  number  and  form  or 
poles  and  zeroes.  As  a  very  optimistic  esti¬ 
mate  of  the  accuracy  of  such  a  procedure,  the 
following  "realistic"  compensation  filter  is 
assumed: 
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The  additional  poles  in  the  filter  are 
required  to  keep  the  gain  at  high  frequency 
finite. 

Using  G  (s)  as  given  by  Equation  (3)  in 
comp 

the  block  diagram  of  Figure  10,  the  pointing 
error  PSD  wns  again  determined  and  is  plotted 
as  Figure  11.  The  RMS  pointing  error  angle 
was  determined  to  be 

0  >13.1  urad 

15  RMS 

It  must  be  emphasized  that  the  above  error 
results  even  though  a  perfect  measurement  of 
the  input  beam  motion,  ,  Is  assumed.  The 
e 


Figure  11. 
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Pointing  Error  PSD  With  Feed  Forward 


error  is  attributable  to  the  imprecision  with 
which  Gcomp(s)  given  by  Equation  (3)  approxi¬ 
mates  the  inverse  of  the  closed  inner  loop 
function  PF^Cs). 

The  base  motion  isolation  performance  of 
the  beam  deflector  elevation  loop  will  be 
further  degraded  for  a  perfect  measurement  of 
input  beam  motion,  ,  is  not  possible.  In 
e 

reality,  such  a  quantity  must  be  measured  by  a 
physical  sensor  uhich  will  add  both  dynamics 
and  noise. 

In  order  to  provide  a  highly  accurate 
measurement  of  ,  a  rate  Integrating  gyro 
e 

will  be  used.  By  providing  a  constant  gain 
feedback  loop  around  such  a  gyro,  the  electronic 
.analog  of  a  mechanical  rate  gyro  is  created 
without  the  mechanical  problems  usually  asso¬ 
ciated  with  rate  gyros. 

The  output  of  a  rate  integrating  gyro  is 
characterized  by 


8 

where 


r  s+l 
g 


h  .  r  .. 

C2,  ¥lUXA“')CMD-VsRA+(v)ugJ'TgU( 


OA 


(4) 


A  «  output  angle  of  the  gyro 
8  I 

Tg  “  -jj8  *  gyro  characteristic  time 

H0  “  8yro  spin  angular  momentum 

>:  ■  gyro  damping 
8 

uIA  >  angular  velocity  along  the  input  axis 
iOqa  “  angular  velocity  along  the  output  axis 
UCMD  "  con,Iianded  angular  velocity 
(v)w  “  gyro  drift  rate  uncertainty  (i.e., noise) 

Since  the  gyro  is  to  be  used  in  a  feedback 

mode,  A  will  remain  very  small.  In  addition, 

*  © 

for  simplicity  of  analysis,  cross-axis  coupling 
twill'  be  neglected  for  the  present.  Thus  the 
rate  integrating  gyro  equation  is  approximately 

H 

_s 

c  . 

Ag  "  [wia-wcmd+<v)ws] 

To  use  the  gyro  as  an  electronic  equivalent 
of  a  conventional  ra^e  gyro,  the  command  angular 
velocity  is  made  proportional  to  A  .  Figure  12 

S 

presents  a  block  diagram  of  such  an  arrangement. 
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Note  that  the  output  is  obtained  by  a  final 
multiplication  by  gain  to  insure  a  unity 

steady  state  gain. 

The  measurement  loop  output,  C(s) ,  is 
related  to  the  desired  input  (s) ,  by 

e 


Ggyro(s)  "  wn  (s) 


■■  £  .8  s2  +  — ja_  o  +  1 

“sKf  HsKf 


The  parameters  for  the  rate  integrating 
gyro  used  to  measure  u0  are  chosen  to  be 
e 

representative  and  correspond  to  those  fot  a 
Kearfott  King  Series  miniature  rate  integrating 
gyro  model  C70  2519  001  4  ,  The  parameter 
values  are 


t  ■  0.006  seconds 
C 

The  choice  of  feedback  gain,  K^,  about 

this  rate  integrating  gyro  i3  of  necessity  a 
compromise.  A  high  gain  produces  a  desired 
increase  in  natural  frequency  but  at  the 
expense  of  j  decrease  in  damping  ratio.  A 
compromise  i  oice  of  «  15  is  made.  This 

choice  result  i  in  a  natural  frequency  of  193.6 
rad/sec  for  the  second  order  dynamics  of  the 
rate  measurement  with  a  damping  ratio  of  0.43. 

Since  the  designar  is  aware  ’that  the  rate 
measurement  process  involves  second  order 
dynamics,  he  can  attempt  to  modify  his  compensa¬ 
tion  filter,  GC0B,p(®)»  to  partially  compensate 

for  the  measurement  dynamics  as  well  as  the 
dynamics  of  the  inner  control  loop.  Again, 
since  the  dynamics  of  the  gyro  measurement 
loop  are  not  known  exactly,  the  feed-forward 


compensation  filter  might  optimistically  be 
chosen  as 


G  (s) 
comp 


n  tzJ_+,1l_a2  +2(0. 7)sx,|  (  s2  j_2(0. 4)Sj_,  I 
°-5(190  1153?  ^iil9P  194~  X) 


(67  +  (5000  + 


(5) 


where  again  additional  poles  have  been  added 
co  provide  finite  high  frequency  gain,  i.e.,  a 
physically  realizable  filter. 


Using  the  above  filter,  the  overall  block  dia¬ 
gram  of  the  elevation  axis  control  loop  is 
shown  in  Figure  13. 

With  G  (s)  given  by  Equation  (5)  the 
comp 

block  diagram  of  Figure  13  was  used  to  determine 
the  pointing  error  PSD  for  the  elevation  axis 
control  loop.  Assuming  no  gyro  noise  in  the 
rate  measurement  loop,  the  RMS  pointing  error 
was  determined  as 

9  -  29.2  urad 

eRMS 

Thus  it  is  apparent  that  the  ability  of 
the  feed-forward  signal  to  compensate  for 
input  beam  motion  is  somewhat  degraded  by  the 
addition  of  measurement  dynamics. 

It  should  be  noted,  however,  that  the 
measurement  dynamics  are  a  function  of  the 
specific  sensor  selected.  Rate  integrating 
gyroscopes  are  available  with  characteristic 
times,  t  ,  an  order  of  magnitude  smaller  than 
8  s 
that  used  in  the  above  calculations  J  .  The 
usefulness  of  such  additional  response  is 
determined  by  the  required  pointing  accuracy 
for  a  specific  mission. 

In  addition  to  dynamics,  an  actual  sensor 
used  to  measure  input  beam  motion,  ,  can  be 

3 
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Figure  13.  Overall  Elevation  Loop 


expected  to  add  measurement  noise.  Adequate 
test  measurements  upon  a  representative  inte¬ 
grating  gyro  to  determine  its  noise  output 
across  the  frequency  spectrum  of  interest  are 
very  scarce.  A  noise  model  PSD  for  a  Honeywell 
gyro  model  CG159C1  has  been  reported  as  6 

Assuming  the  gyro  noise  to  be  uncorrelaced 
with  all  other  disturbances,  the  performance  of 
,  the  elevation  control  loop  was  again  evaluated 
using  the  above  gyro  noise  PSD.  The  RMS  point¬ 
ing  error  caused  by  gyro  noise  alone  was  deter¬ 
mined  to  be 

0  (syro  no'ise)  »  0.19  ur&d 
CRMS 

Recently,  information  on  the  noise  PSD's 
of  several  high  quality  gyroscopes  has  been 
published  The  noise  is  presented  in  terms 
of  equivalent  input  angle  in  arcsech  Across 

hi; 

a  spectrum  of  0.01  hz  to  100.00  hz,  several 
of  the  tested  gyros  are  shown  to  possess  RMS 
uncertainties  of  less  than  0.03  arcsec  of 
lnpuc  angle.  This  is  consistent  with  the  RMS 
pointing  error  determined  above  and  indicates 
that  short  term  gyro  noise  is  an  unlikely 
source  of  significant  error  in  the  operation 
of  precision  optical  directors. 

V.  AZIMUTH  LOOP  DESIGN 

The  basic  azimuth  loop  of  Figure  5b  is 
redrawn  as  Figure  14.  It  is  seen  that  the 


rotation  of  the  pointing  direction  about  the 
p  axis  results  not  only  from  rotations  about 
the  k  axis  which  are  controlled  but  also  from 
rotations  about  the  n  axis.  Inasmuch  as  the 
controlled  direction  has  no  gimbal  Isolation 
from  the  base  about  the  n  axis,  rotations  of 
the  base  about  this  axis  produce  rotational 
disturbances  of  the  line-of -sight. 

In  a  manner  analogous  to  that  described 
in  the  preceding  sections,  a  feed-forward  path 
can  be  added  to  reduce  the  effect  of  base  motion. 
Such  a  loop  is  shown  in  Figure  15.  The  angular 
acceleration  of  the  base  about  the  n  axis  is 
measured  and  provided  as  a  command  torque  to  the 
azimuth  torque  motor. 

The  performance  function  for  the  angular 
accelerometer  corresponds  to  that  measured  for 
a  Systron-Donner  Model  4590-F-l-AC  angular 
accelerometer  ?  . 

As  with  the  elevation  control  loop,  a 
filter  is  added  to  compensate  for  the  measured 
dynamics  of  the  control  loop  and  of  the  angular 
accelerometer.  Again,  the  matching  of  the 
filter  to  the  actual  dynamics  is  chosen  to  be 
on  optimistic  estimate  of  that  which  could  be 
attained. 

Using  the  base  motion  PSD  of  Equation  (2) 
and  the  above  control  loop,  the  PSD  of  azimuth 
pointing  error  about  the  p  axis  is  determined. 
The  RMS  value  of  this  pointing  error  produced 
by  base  motion  is  found  to  be 

9  ■  14.2  ur.id 

p  RMS 
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Figure  1A.  Basic  Azimuth  Loop 


Figure  15.  Azimuth  Loop  With  Feed  Forward 


VI.  ADDITIONAL  CONSIDERATIONS 

In  an  attempt  to. assess  the  minimum  point¬ 
ing  error  which  might  be  attainable  with  an 
airborne  two-axis, beam  deflector  of  the  design 
considered,  several  complications  have  been 
neglected.  These  include  gyro  output  axis 
cross-coupling,  gimbal  bearing  friction,  and 


gimbal  non-rigidity. 

As  indicated  by  Equation (A), an  angular 
velocity  of  the  gyro  case  about  its  output  axis 
will  produce  an  output  signal  from  the  gyro 
which  is  not  desired.  In  the  control  loop  pre¬ 
sented  for  elevation  axis  control,  it  is  assumed 
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that  a  single  degree  of  freedom,  rate  Integrating 
gyro  is  affixed  to  the  beam  deflector  with  its. 
input  axis  aligned  along  the  e  axis.  The  azimuth 
control  loop  assumes  a  similar  gyro  mounted  on 
the  outer  gimbal  with  its  input  axis  along  the 
k  direction.  In  order  to  avoid  output  axis 
coupling  between  the  loops,  the  output  axis  of 
this  azimuth  gyro  can  be  aligned  with  the  n 
direction.  On  the  other  hand,  the  output  axis 
of  the  elevation  gvro  must  lie  in  a  plane  con¬ 
taining  the  k  direction.  Hence  controlled 
angular  velocities  of  the  azimuth  loop  about 
the  k  axis  will  produce  an  erroneous  output  of 
the  elevation  gyro  due  to  its  output  axis  cou¬ 
pling.  This  coupling  is  diagrammed  in  figure 
16.  Although  such  coupling  provides  an  unwanted 
disturbance,  it  does  not  alter  the  loop  stability. 


Figure  16.  Gyro  Cross  Coupling  Between  Loops 


In  addition  angular  velocities  of  the  base 
about  the  n  direction  provide  erroneous  outputs 
from  both  the  azimuth  and  elevation  loop  gyros 
due  to  output  axis  coupling. 

The  effect  of  output  axis  coupling  can  be 
reduced  by  compensation.  For  example,  recog¬ 
nizing  that  angular  velocities  about  the  n  and 
k  axe3  will  couple  to  the  elevation  control 
loop  gyro,  one  can  measure  these  angular  veloc¬ 
ities  with  appropriate  gyros.  Knowing  the 
characteristic  time  of  the  elevation  gyro,  a 
compensating  command  torque  can  be. provided 
to  compensate  for  the  unwanted  output  axis 
coupling. 

In  addition  to  the  effect  of  gyro  cro3s-axis 
coupling,  the  problems  of  gimbal  bearing  friction 
and  gimbal  non-rigidity  have  not  been  included 
in  this  analysis.  These  effects  could  alter 


the  natural  frequencies  and  stability  of  the 
loops.  Each  of  these  neglected  effects  can  be 
expected  to  increase  the  pointing  errors.  The 
amount  of  degradation  of  loop  performance 
could  depend  upon  the  specific  mechanical 
design  of  both  the  beam  deflector  an'  its 
gimbals. 


VII.  RESULTS 

For  the  control  loop  design  end  base  motion 
spectrum  considered  in  this  investigation,  the 
RMS  pointing  error  for  the  cases  considered  is 
summarized  as  follows: 

ELEVATION  LOOP 

8  in  prad 

-HMS _ 

1.  No  feed-forward  compensation  9A2.8 

2.  Feed-forward  compensation*. 

Perfect  mcas\:reme>'.r.  of  13.1 

e 

3.  Feed-forward  compensation*. 

Gyro  2nd  order  dynamics  in 
measuring  ioQ  **.  29.2 

e 

A.  Feed-forward  compensation. 

Gyro  2nd  order  dynamics*  and 

$yro  noise  in  measuring  w  29.2 


AZIMUTH  LOOP 

9  in  prad 
PRMS _ 

Feed  forward  compensation.* 

Accelerometer  2nd  order  dynamics 
in  measuring  .  1A.2 

n 


VIII.  CONCLUSIONS 

1.  The  use  of  feed-forward  compensation 
can  significantly  reduce  the  pointing  error 
caused  by  base  motion  of  an  airborne  two-axis 
beam  deflector. 

2,  The  resulting  pointing  error  depends 
upon  the  degree  with  which  hardware  dynamics 
can  be  measured  and  matched. 


*  Feed  forward  compensation  includes 
"optimistic"  matching  of  control  loop  dynamics. 

**  Effect  of  gyro  dynamics  can  be  induced  by 
choosing  a  gyro  with  a  smaller  characteristic 
time. 


3.  Assuming.- the  use  of  a  high  quality 
gyro,  noise  in  the  gyro  used  to  measure  base 
motion  for  the  feed-forward  compensation  does 
not  significantly  increase  pointing  errors. 


Derivation  of  Dynamic  Equations 

This  section  is  intended  as  a  detailed 
derivation  of  the  dynamic  equations  of  the 
two-axis-girl) ailed  beam  deflector  as  depicted 
ir.  Figure  2. 

The  following  definitions  are  used  in 
the  development: 


axes  and  that  the  gimbals  are  balanced  about 
their  respective  axes  of  rotation. 

Inner  Glmbal  Dynamics 


INERTIAL 


E -Torques  on  the 

inner  glmbal  (A.l) 
about  its  c.m. 


Applying  the  law  of  Coriolis 

dt  “  dt  +  V^JmV 

INERTIAL  M 

a 

■JM  d?|  +  V(ImV 


Coordinate  Frames 

V  frame:  i,j,k  coordinates  fixed  to  the 
aircraft 

0  frame:  n,e,k  coordinates  fixed  to  the 
outer  glmbal 

M  frame:  r,e,d  coordinates  fixed  to  the 

mirror,  i.e.,  to  the  inner  glmbal 


Writing  coordinates  in  the  I  frame 


Af  0  0 


0  A  0 

e 

0  0  Aj 


INERTIAL  frame:  inertial  frame  with  unspecified 
coordinates  axes 

Notation 


"  « 

M  M_ 


Superscript  indicates  the  coordinate  frame  in 
which  a  vector  is  coordinated. 


Subscript  Indicates  the  coordinate  frame  with 
regard  to  which  derivatives  are  considered. 

Varlebles 

n  ■  Mirror  azimuth  angle 


dt  UM  "  “M_ 


e  ■  Mirror  elevation  angle 


Uy  ,  U  ; 

i  j 


“o  ’  “o 


wv  "  Aircraft  angular  rates 
vk  about  i,j,k  axes, 
respectively. 

»  Outer  glmbal  angular 
uk  rates  about  n,e,k  axes, 
respectively. 


«M  ,  *  Inner  glmbal  (mirror) 

r  e  "d  angular  rates  about  r,e,d 
axes,  respectively. 

B^,  B  ,  B.  ■  Principal  moments  of  inertia 
n  e  for  outer  glmbal  about  n,e,k 

axes,  respectively. 

A  ,  A  ,  A ,  ■  Principal  moments  of  inertia 
e  for  inner  glmbal  (mirror) 

about  r,e,d  axes,  respectively. 


vVhI  ■ 


UM  UM  V 

r  <*,  a 

ArwM  Ae“M  Ad“M, 

red 


(Ad"Ae)  \  \ 
d  e 

<VAd>  \  wMd 

(VV  \ 


It  is  assumed  that  glmbal  axes  are  principal 


Substituting  (A. 3)  and  (A. 4)  into  (A.2) 

fAr  \  +  ‘W  UMd  \ 


Torques  -  A  (^  +  (Ar-Ad)  uM  (A.  5) 

*  e  r  d 

A,  V  +  (Ae-V  “M  “M 
_  a  a  e  rJ 


BOO 

n 

-  0  0 
o  e 

0  0  B, 


Taking  the- component  of  (A. 5)  along  the  e  axis 
yields  (neglecting  bearing  friction) 

Te  "  Ae  “H  +  (Ar'V  %  “m  (A'6) 

e  dr 

however,  due  to  gimbnl  constraint: 


uo’'  “o 


UM  ' 

'  w0 

cos  c  -  <0.  sin  e 

(A.  7) 

r 

n 

°k 

d 

-  1° 

“h, 

’  W0 

sin  e  +  cos  e 

(A.  8) 

dt 

«0 

d 

n 

°k 

0 

UM  UM .  " 

<*0 

-  u2  )  sin  e  cos  e 

r  d 

n 

°k 

+  wo 

(cos2  c  -  sin2  e) 

(A. 9) 

n 

n  k 

0 

Substituting  (A. 9) 

into  (A. 6)  yields 

“o  x  Vo  “ 

u0 

n 

AeuM  "  Te 
0 

-(A, 

rV  (wo  i sln  2 

n  k  2 

C 

B  u. 
n0n 

-  (A  -A.)  u  uQ  cos  2  r  (A. 10) 

it  k 

equation  (A. 10)  describes  the  dynamics  of  the 
inner  global, about  the  elevation  axis  (e  axis) 
as  a  function  of  outer  global  angular 
velocities. 

Outer  Global  Dynamics 


INERTIAL 


_  —^Torques  on  the 

(J-U  )  -  }  outer  gimbal  (A. 11) 
u  0  “about  its  c.n. 


The  torques  on  the  outer  gimbal  are  pro¬ 
vided  by  the  azimuth  torque  motor,  reaction 
from  the  elevation  torque  motor,  and  by  gimbal 
restraint  from  the  base  and  from  the  inner 
glmoal. 

£  <Vo>  ■  T0  It  “o  +  “o  x  Vo 

INERTIAL  0 

(A.  12 

Writing  coordinates  in  the  G  frame 


<YV  uo  \ 

e  k 


“0  X  Vo  '  <W  u0,  u0  (A>16) 

k  n 


(Be-Bn)  o0  m0 
n  e 


Substituting  (A.  15)  and  (A. 16)  into  (A. 12) 


Bn  “0  +  (Bk'Be)  w0  “0. 
n  e  k 


^Torques  -  Bfi  +  (B^-B^;  W()  U(.  (A.  17) 
o  k  n 


Bk  X  +  <W  w0  “o 

k  n  e 


Consider  the  torques  on  the  outer  gimbal 
along  the  k  axis.  These  are: 

a)  The  k  component  of  the  torque  T,,. 

MO 

exerted  by  the  inner  gimbal  on  the  outer 
gimbal. 


b)  T^»  the  azimuth  control  torque.  Now 
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t  -  Jr  where  T  ‘  is  the  torque  exerted  by 
lM0  CM  OM  ~  . 

the  outer  global  on  the  Inner  global.  T0M  is 

given  ^y  Equation  (A.5)  In  H  coordinates. 

Thus  the  component  of  along  the  V. 
axis  is : 

*»  •  V  -Ad%  008  6  +  (  W  \\  C0S  C 

d  (A.  18) 

,  +AruH  sin  c  +  (Ad"Ae)  «M  “«  sin  e 


-  d 

-  ,  d 

“M  "  dt 

“0+  d? 

0 

0 

Writing  (A.2S)  in-M  frame  components 

0 


—  “  "  4:  “n  +  “m  ”  “o  (A.26) 

dt  M  dt1  0  “e  °e 

0  0 


Substituting  (A. 18)  into  (A.17) 

Vok“(W\\  +  Tk  (A' 

+  Sin  c  +  (Ad-Ae)»Me«Md  sln  E 

‘  Vm  C0S  6  +  (Ar'Ae)uMeUMr  C°S  6 


Substituting  (A.22)  into  (A.26)  yields 


d 

dt 

WM 

0 

u„  cos  e  -  to,  sin  e 
0 

n  * 


sin  e  +  Wq  cos  e 


&  %  "  \ 
M 


d  - 

— _  Q  (4) 

dt  o  0 


From  (A. 21) 


dFl  w0 


w„  cos  e  -  wn  sin  e 

U  v't» 

n  * 


w  sin  e  +  w0  cos  c 


, - M  -M 

<vV  XUM 


“o  ”WM 
e  e 


«0  cost  -  to  Sint  *  to^sinc  +  topcoat) 
n  " 

Substituting  (A. 28)  and  (A. 27)  into  (A.-3) 
yields 

i.  “  to0  cost  -  to0  Sine  +  toQ  toQ  sine  +  to  cost 
nr  un  k  e  n  e  * 

-  sine  -  «iM  to  cost  (A.29) 

e  n  c  k 


£  V  x “m 


Similarly 

“u  "  “n  sinE  +“o  C0SE  "  W0  “0  c08e 
Md  n  k  e  n 


+  u„  wQ  sint  +  uM  toQ  cost  -wH  “0  sir 
'‘e  k  e  n  e  k 

Substituting  (A.29)  and  (A. 30)  into  (A. 19) 
yields 


(A.  30) 
sine 


How,  obviously 


WM  "  “o  +  (“M  ‘  “0) 
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\  “ot. 


(B  -B  )  u.  un  +  T. 
n  e  0  0  k 
n  e 

+A  i.  cose  sine  -  An.  sin2'  +  A  <on  sin2e 
r  o  r  u,  r  u  u 

n  k  e  n 

+A  <i>0  uQ  sine  cose-A  u.,  sin2e-A  uq  sine  cose 
r  e  k  r  e  n  e  k 

-  A^uq  sine  cosc  -  c'.)s2e  +  A^Ug  uq  cos2e 

n  e  n 


Adu0  uq  sine  cose  -  Aji*^  Uq  cos 2e 
e  k  e  n 


+  Ad<uM  uQ  sine  cose  +  AjU^  sine 
e  k  e  d 


+  A  u,,  cose  -  A  u„  (w„  sine  +  u„  cose) 
r  n  m  e  M  M ,  ii 

e  r  e  d  r 


Now 


<i)u  ■  u_  sin  e  +  (i).  cos  e 
M.  U  U. 

d  n  k 


(A.31) 

(A.  32) 
(A. 33) 


ul,  “  uQ  cos  c  -  “0  sin  e 
r  n  °k 

Substituting  (A.32)  and  (A. 33)  into  (A.31) 
yields 

(Bk+Arsir.2e  +  Adcos2e) 

k  'i 

+  (Bn-Be)u0  aQ  +Aru)Q  cose  sine  +  A^Ug  u0  sin2e 
n  n  e  n  e  n 


Again  rearranging  (A. 34)  produces 

(Bj+Arsin2e  +  Adcos2e)uQ  » 

+  (B  -B  +A  8in2e  +  A.cos2e)u.  u- 
ner  d  '0  0. 

n  e 


-  (A  -A  cos  2  e  +  A. cos  2  c)u.  w.. 
e  r  a  On 

n 


+  (A*-A  )  'sin  e  cos  e  2  «tl 

dr  «e  °k 

-  (Ad'Ar)  sin  e  cos  e  uQ  uQ 

-  (A.-A  )  sin  e  cos  e  u. 

a  r  u 


Now, 


L\ 


d 

V 

i 

• 

dt 

“v 

m 

V 

V 

i 

“V. 


+  A_w„  «A  sin  e  cos  c  -  A  uu  u.  sin2e 
too,  r  n  u 

e  k 


c  n  u  \ 

.  .  \ 


-  A  iij  sin  e  cos  e  -  A^Wg  sin  e  car.  e 

r  e  k  n  9 

+  A ,w„  un  cos2e  -  A,u_  un  sin  e  cos  e 
a  0  u  a  o  o. 

e  n  e  k 

-  A.Wu  un  cos2e  +  A.w,.  un  sin  e  cos  e 

a  M  o  and. 

e  n  e  k 

+  A,wu  (u„  sin2e  +  o>.  sin  e  cos  e) 
a  M  o  o, 

e  n  k 

+  A  uu  (u.  co82e  -  u>.  sir.  e  cos  c) 
r  pi  o  u. 

e  n  k 


"  Ae  WH  “0 
e  n 


\ 

\ 

\ 

<A\34) 

\ 

\ 


\ 


d 

dt 

«V 

V 

p  uv  cos  n  +  sin  n" 

i  j 


sin  n  +  cos  n 
i  3 


V 


Now 


0  dt 


u0  +  (wv”“o)  X  “o 


Again,  obviously 


“0  ’  “V  +  (w0  “  UV) 


Thus 


_d_ 

dt 


d 

~  .  ± 

w0  '  dt 

WV+  di 

V 

V 

(VV 


■“V 


(A. 35) 


(A. 36) 


(A. 37) 


(A.  38) 


(A. 39) 


(A. 40) 


(A. 41) 


Writing  (A. 41)  in  0  frame  components 

.  .  .  To 


d 

0 

_d_ 

dt 

“0 

m  , 

dt 

“v 

V 

V 

.  [V  \j 

Substituting  (A. 38)  into  (A. 42)  yields 


d  - 
dt  “o 


w„  cos  n  +  m„  sin  n 

i  J 


sin  n  +  m„  cos  n  (A.  43) 
vi  vj 


Equation  (A. 47)  describes  the  dynamics  of  the 
beam  deflector  about  the  k  axis  in  terms  of 
base  motion  and  motion  of  the  beam  deflector 
about  the  k  and  e  axes. 


Appendix  B 

Disturbance  Input  Relationships 

The  disturbance  inputs  to  the  elevation 
loop,  T.  and  w.  ,  and  those  to  the  azimuth 
e  e 

loop,  and  w.  ,  can  be  written  in  terms  of 
a  °n 

base  motion  and  required  line-of-sight  rates 
by  noting! 

w.  "  cos  n  +  w  sin  n  (B.l) 

°n  vi  vj 


_  _  o  _o 

<«V  “  V  x  “  0 


(A. 44) 

k  I 


v\ 

k  k 


cosn+Wy  sinn  sinn+M  coan  u.. 

vi  vj  vi  vj  Vk 

Substituting  (A. 44)  and  (A. 43)  into  (A. 40) 
produces 

•  •  • 

wo  “  wv  cos  n  +  wv  sin  n 

n  i  J 

+  Uy  uv  sin  n  -  w  My  cos  n 
vi  k  vj  vk 

-  sin  n  +  My  u0  cos  n 

vi  \  vj  uk 


w.  ■  -  My  sin  n  +  u„  cos  n  (8.2) 
e  vi  vj 

Substitution  of  (B.l)  and  (B.2)  into 
(A. 10)  and  (A. 47)  results  in  all  disturbances 
expressed  in  term3  of  base  motion  rates  and 
required  line-of-sight  tracking  rates. 

In  order  to  determine  the  I.'SD  of  the 
pointing  error,  it  is  necessary  to  know  the 
PSD's  of  the  various  disturbance  Inputs.  This 
can  be  readily  accomplished  knowing  the  PSD  of 
each  component  of  base  motion  and  required 
tracking  rates  if  Independence  is  assumed 
amongst  the  various  angular  velocity  components 

As  an  example  of  the  determination  of  PSD 
consider  the  disturbance  input  term  u  au  . 

Vi  Vj 

The  PSD  for  My  or  My  is  give.j  by 

i  J 

Equation  (2). 


a.  “  -  m  sin  n  +  My  cos  n  (A. 46) 
e  vi  vj 

Substituting  (A. 45)  and  (A. 46)  into  (A. 35) 
yields 

(B.  +  A  sin2e  +  A,  cos2e)  w.  ■ 
k  r  d  0k 

Tfc  +  (Bn-Be+Ar  sin2c  +  A.  co32e)  uQ  uQ 
t.  -  n  ye 

-  (A -A  cos2e  +  A.cos2e)m-  mw  +(A.-A  )sint  cose  2m„  m„ 

e  r  d  0nMe  dr 

-  sin  c  cos  e  jwv  cos  n  +  My  sin  n  +  «v  My  sinn 


-  My  w„  cos  n  -  2m„  u-  sin  n  +  2m, ,  w.  cos  n 
Vj  Vk  vi  °k  M  °k 


Assuming  independence,  it  is  directly 
shown  that 

R  (t)  -  R  (t)  R  (t)  (B.3) 

“v  wv  “v  wv 

vi  J  vi  vj 

where  R  (t)  is  the  autocorrelation  of  the 
z 


variable  z. 

Taking  the  inverse  Fourier  transform  of 
Equation  (2)  yields 

(  -1.5|xj  -1400|t| 

\  (T)  “  \  (T)  *  0,1  \  3  ^7800 - 


Substituting  (B.4)  into  (B.3)  and  taking 
the  Fourier  transform  yields 
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